A Quarkyonic Phase in Dense Two Color Matter? 
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We present results from simulations of Two Color QCD with two Wilson quark flavors in the 
presence of a quark chemical potential jj, at two different lattice spacings. The equation of state, 
conformal anomaly, superfluid order parameter and Polyakov line are all discussed. Our results 
suggest that the transition from hadronic to quark matter, and that from confined to deconfined 
matter occur at distinct values of /x, consistent with the existence of a quarkyonic phase in this 
model. 
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A quantitative understanding of cold dense baryonic 
matter, ie. the region of the (T, fi) plane convention- 
ally placed at the lovifer right of the QCD phase diagram 
(here /x is the quark chemical potential) , remains an out- 
standing problem in theoretical physics. There has been 
much discussion of possible exotic color superconduct- 
ing (CSC) phases in this region, where color-carrying de- 
grees of freedom such as quarks or gluons are all gapped 
via a Higgs-Meissner mechanism, and which may also 
be superfluid or even crystalline This picture has 
its firmest theoretical support at weak gauge coupling, 
where superconductivity is understood as being due to 
the BCS condensation of diquark Cooper pairs. More 
recently, an alternative scenario based on large-A^c ar- 
guments has emerged, whose ground state at sufficiently 
high density has a restored chiral symmetry, so that ther- 
modynamic properties are well-described by a degenerate 
relativistic Fermi sea of quarks characterised by a Fermi 
momentum kp, but which remains in a confined "quarky- 
onic" phase so that all excitations necessarily carry color- 
singlet quantum numbers 

It remains unclear whether the CSC and quarkyonic 
pictures are truly distinct, or are in some sense comple- 
mentary. However, it is salutary to recall that even basic 
questions involving dense matter, such as the maximum 
stable mass of a neutron star, need quantitative input 
about the equation of state of ultradense matter regard- 
less of the nature of the ground state; this requires a 
controlled non-perturbative calculation. 

The most reliable source of such information, lattice 
QCD, is in general inoperable in this regime for the fol- 
lowing reason. In Euclidean metric the QCD Lagrangian 
density for quarks reads 
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It is straightforward to show 75M(/i)75 = M''{—fi), im- 
plying detAjf(/i) — (detM(— /x))*, and therefore that the 



path integral measure is not real and positive for n 0. 
Monte Carlo importance sampling, the mainstay of nu- 
merical lattice QCD, is ineffective. It is helpful to con- 
sider what goes wrong if the real positive measure factor 
deiM'^'M as implemented in, eg. the hybrid Monte Carlo 
(HMC) algorithm, is used. In QCD, while M describes 
a color triplet of quark fields q, the factor describes 
color antitriplet "conjugate quarks" q'^. Gauge singlet 
bound states of the form qq'^ resemble mesons, but carry 
non-zero baryon charge B > 0. The lightest such state 
is degenerate with the pseudo-Goldstone 7r-meson; hence 
HMC simulations with fi predict an unphysical "on- 
set" transition from the vacuum to a state with quark 
density riq > at fio — ^itLt^. The resulting ground state 
is a Bose-Einstein condensate (BEC) of diquark baryons, 
bearing no resemblance to nuclear matter, which phe- 
nomenologically we know forms at /io ~ ■^mnucieon- 
The physical transition can only be found if the correct 
complex path integral measure det^Af is used, and must 
result from extremely non-trivial cancellations between 
configurations with differing phases - this has come to 
be known as the Silver Blaze problem 3i|. 

In this paper we consider an alternative strongly- 
interacting theory. Two Color QCD (QC2D), in which 
the gauge group is SU(2). Since q and q live in equiva- 
lent representations of SU(2), it follows that detM(/i) = 
detr2M*(/i)T2 is real and therefore the theory has a 
positive measure for an even number Nf of quark fia- 
vors 0]. Physically this is expressed through both qq 
mesons and qq, qq baryons falling in the same hadron 
mulitplets. For sufficiently light quarks the scale hierar- 
chy <Si mp permits the use of chiral perturbation the- 
ory (xPT) in studying the response of the lightest multi- 
plet to ^ ^ [HI. The key result is that for /i > /Xq = ^m^ 
a non-zero baryon charge density > does develop, 
along with a gauge-invariant superfluid order parameter 
which for Nf = 2 reads (qq) ^ (V'*''C'75'^2eab^) 7^ 0, 
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where T2 acts on color indices and Cab = —^ba on flavor. 
The resulting BEC is composed of weakly interacting qq 
baryons with = 0+. 

For > Ho leading-order xPT predicts a smooth ro- 
tation of the chiral condensate (V'^) into the superfluid 
condensate {qq) as /i increases [5[- In addition there is a 
quantitative prediction for quark density: 



(2) 



where the parameters /io and suffice to specify xPT at 
this order. It is possible to develop the thermodynamics 
of the system at T = more fully, to extract pressure 
and energy density [6]: 



J u„ 



AN,fl[^?^i^^-2^,l]■ (3) 



Note that the trace of the stress-energy tensor 
iTf^fj^)^PT < for fi > y/3^io- 

These model results should be contrasted with those 
of another paradigm for cold dense matter, namely a de- 
generate system of weakly interacting (thus presumably 
deconfined) quarks populating a Fermi sphere up to some 
maximum momentum kp fa Ep = fi: 



£SB = SPSB 



(6) 



In this picture superfluidity arises from condensation of 
diquark Cooper pairs from within a layer of thickness A 
centred on the Fermi surface; hence 



(qq) cx A/i^. 



(7) 



Since psp eventually exceeds PxPT as /x increases, the 
degenerate system must be the more thermodynamically 
stable at high density. In fact, since both and e are 
discontinuous at the point where pss — PxPT, this naive 
treatment predicts the resulting deconfining transition is 
first order [6[. 

These considerations have motivated us to pursue lat- 
tice simulations of QC2D beyond the BEC regime, using 
Nf = 2 flavors of Wilson fermion. The quark action is 

S=Y^ i^.M^^ + KjlVr (^75)^2^1 - h.c], (8) 



i=l,2 



with 



+ (l + 7.)e"^'"'C/^(y)'^y,.-^l. (9) 



A conventional Wilson action was used for the glue fields. 
Further details can be found in 

Since Wilson fermions do not have a manifest chi- 
ral symmetry, we have little to say at this stage about 
this aspect of the physics, which at high quark density 
should be of secondary importance for phenomena near 
the Fermi surface; they do however carry a conserved 
baryon charge due to the U(1)b symmetry ^ e^^Vj 
-0 I— ■06"'". Our initial runs on a 8'^ x 16 lattice with 
/? = 1.7, K = 0.178 corresponding to lattice spacing 
a = 0.230(5)fm, m^a = 0.79(1) and m^/nip = 0.779(4) 
have been described in . In this Letter we present data 
from runs on an approximately matched 12^^ x 24 lattice 
with /3 — 1.9, K = 0.168 corresponding to a = 0.186(8)fm, 
m-^a — 0.68(1) andm^/mp — 0.80(1). The physical scale 
is set by equating the observed string tension at ^ = 
to (440MeV)^. Note that the physical temperature T 
is approximately 54(l)MeV for the smaller lattice, and 
44(2)MeV for the larger. We used a standard KMC al- 
gorithm - the only novelty is the inclusion of a diquark 
source term (proportional to j in Eqn.®); this mitigates 
the impact of IR fluctuations in the superfluid regime and 
also enables the algorithm to change the sign of detM for 
a single flavor, thus maintaining ergodicity. All results 
presented here were obtained with ja — 0.04; ultimately 
the physical limit j — >■ must be taken. 
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FIG. 1: (Color online) rig/nsB and p/psB vs. fj, for QC2D. 
Inset shows Sq/esB for comparison. 

Fig. [T] shows results for quark density and pressure 
as functions of /i, plotted as ratios of the same quan- 
tities evaluated for free massless quarks on the same lat- 
tice 0. In the J — >■ limit the onset is expected at 
/iofl = 0.34 corresponding to fio — 360MeV; the obser- 
vation oi nq,p > for fi < fio is an artifact of working 
with j ^ 0. Beyond onset the ratio riq/nsB rises to a 
peak at /i « 400MeV, then falls to a plateau beginning 
at fiQ « 530MeV, which continues until hd ^ 850MeV 
where it starts to rise again. If following the arguments 
presented above we associate the plateau with the setting 
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FIG. 2: (Color online) Conformal anomaly T^f^a"^ vs. jj., plot- 
ted together with separate quark and gluon contributions. 



in of degenerate matter then we identify a BEC/BCS 
crossover at ~ 530MeV, corresponding to a quark 
density « 4 — 5fm~^, ie. roughly 10 times nuclear 
density. In order to specify these numbers with greater 
precision, as well as the statistical errors manifest in the 
error bars of Fig. [T]we need to establish control over lat- 
tice artifacts by taking the continuum limit a — ^ 0. So 
far we have data from only two lattice spacings; results 
for Uq/nsB taken with a — 0.23fm plotted in Fig.[l]show 
that there is reasonable scaling for /i ^ 400MeV - in- 
deed, the small difference in the physical j between the 
two ensembles has maximal impact for /i « @, ■ The 
situation at larger /i will be discussed further below. 

In contrast to xPT, the quark contribution to the en- 
ergy density Sq exceeds the free field value by almost a 
factor of 20 for /i ^ /Xq, as shown in the inset of Fig.[TJ it 
should be remarked here that unlike Uq and p, e is sub- 
ject to a multiplicative quantum correction known as a 
Karsch coefficient 12] which is still to be calculated for 
this system, though its renormalised value is unlikely to 
differ by more than 50%. In any case, since the Karsch 
coefficient is /i- independent, the shape of the curve will 
remain the same. Because of this unexpected behaviour 
at small the energy per quark Sq/uq exhibits a shallow 
but robust minimum for /i > ^q, a feature completely 
absent in the model governed by Eqns. ([2][6]). 

Fig. [2] plots the conformal anomaly T^^ — {T^^)g + 
{Tf,f,)q, with 
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(10) 
W)), (11) 



LCP 



where LCP denotes the beta-function is determined 
along a line of constant physics, and once again a vac- 
uum contribution must be subtracted. The required 
beta-functions are estimated from the parameters of the 



two matched lattices to be -0.85(17) (.g) and 0.042(9) (g) 
(these errors are not included in Fig. [2])- 

Fig. [5] shows that for /i ^ /ig the contribution of both 
quarks and gluons to T^^ is initially positive and in- 
creasing with /i; here both quarks and gluons are con- 
tained within tightly-bound non-relativistic bosons, im- 
plying that e > ?>p. For /i ^ fig, however, their be- 
haviour diverges sharply (see inset), which could possibly 
be explained by quark and gluon degrees of freedom now 
being governed by differing quantum statistics. In fact, 
the gluon data are very well approximated over the whole 
/i-range by a parabola, and (T^;^)g accordingly becomes 
negative for /i ^ /ijj. This change in sign has also been 
predicted using xPT and the property of asymptotic free- 
dom [7] (another way of understanding the necessity for 
{Tfj.fi)g to change sign is that the plaquette must revert 
to its quenched value in the limit /i — >■ oo where quantum 
corrections due to quarks are Pauli-blocked fej). At the 
same point there is a very sharp change in the behaviour 
of which had been approximately constant for 

^ fJ. ^ fJ,D- Fig. [2] shows that at large fi the quark 
contribution dominates, so that lim^_j.oo > 



behaviour is not predicted by xPT (eg. Eqn 
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Ref. 0), although the positivity of T^^ in this limit is 
consistent with three- loop perturbation theory ■ 
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FIG. 3: (Color online) Gluon energy density Eg/fi'^ versus fi. 



Next consider the gluonic energy density given by 
3/3 



(12) 



(once again, a /j,-independent Karsch coefficient is still 
to be determined). Fig. [3] plots the dimensionless com- 
bination Eg/ jjL^ against /i; of course this quantity is not 
predicted either in xPT or the free quark gas. While we 
therefore have no quantitative theory of the gluonic con- 
tribution to QC2D thermodynamics at /i 7^ 0, we would 
expect its relative importance to increase across a decon- 
fining transition. In fact, the ratio is remarkably constant 
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over a wide range of fj,, consistent with dimensional anal- 
ysis; in particular there is no sign of singular behaviour at 
/i = fiQ, although there is a systematic rise for /ijj. 
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FIG. 4: (Color online) Superfluid order parameter {qq) / 
and Polyakov line versus ^. 

Fig.|4]plots quantities giving information on the nature 
and symmetries of the ground state. In the limit j — > 0, 
the diquark condensate {qq) is an order parameter for the 
spontaneous breaking of U(1)b leading to baryon num- 
ber superfluidity. Although the data of Fig. |4] are taken 
with j ^ 0, implying some care must be taken with the 
extrapolation j — >■ at small /i [W| , we are confident that 
this symmetry is broken for all /i > /io. The approximate 
flatness of the curve for /xg ^ ^ /i/j is then evidence 
for a scaling {qq) cx /i^ similar to Eqn.([7]). We take this 
as an indication that in this region the system consists of 
degenerate quark matter with a Fermi surface disrupted 
by a BCS instability. 

The Polyakov line is an order parameter for decon- 
finement in the limit of infinitely massive quarks ~ away 
from this limit it continues to yield information on the 
free energy of an isolated color source. Fig. |4] shows that 
QC2D remains confined for fi < fiz)^ but that there ap- 
pears to be a transition to a deconfined state for chem- 
ical potentials in excess of this value. In physical units 
fj-D ~ 850MeV, corresponding to quark density nq « 16 
- 32fm^^, some 35 - 70 times nuclear density. 

To summarise, the simulations suggest that QC2D has 
three distinct transitions (or at least crossovers). The 
first, at /i = /io, is a firmly established second order 
phase transition (in the limit j — )• 0) from vacuum to a 
BEC superfluid, and is described accurately for the most 
part by xPT (the quark energy density Sq looks to be 
an important exception) . Since the pion is not especially 
light with our choice of lattice parameters, implying only 
a moderate separation of Goldstone and hadronic mass 
scales, the /^-window within which the BEC is favoured 
is not particularly wide. 

The second transition at /x = iiq looks like a BEC/BCS 



crossover to form a ground state where the scalings of the 
observables ng(/i), p(/i), £qi^J■) and {qq{^)) all suggest it 
is formed of degenerate quark matter with a well-defined 
Fermi sphere, albeit one whose surface is disrupted by 
a BCS condensate. We note that effective treatments 
based on both xPT and NJL models predict Uq/nsB to 
be monotonic decreasing in this regime, and are unable 
to fit the lattice data [9|. The distinct nature of this 
region is also supported by the diverging behaviours of 
{Tf_,f_,)g and (r^^)g for n Z hq seen in Fig. [21 ahhough 
the reason for the peculiar behaviour of {T^^)q is not well 
understood at present. The transition at /i ~ /xq is most 
likely a smooth crossover, but the exact nature of this 
putative transition requires further study, in particular a 
careful extrapolation to the limit of zero diquark source. 

The third transition at /i = /iu is signalled by a change 
in the scaling of the thermodynamic observables, notably 
£g(/i) and {T^^)q, a change in the sign of {T^^)g, and a 
non-zero Polyakov loop. For /i > /i^ the system consists 
of deconfined quark matter. 

An immediate concern is the validity of the deconfin- 
ing transition in the continuum limit. Fig. |4] shows that 
with a = 0.23fm ~ 600MeV, and is practically indis- 
tinguishable from HQ ; for this reason only a deconfined 
quark matter phase was identified in Ref. [6] (Cf. Fig.[T]). 
In both cases, however, the quark density in lattice units 
Uqa^ = 0.17 (coarse) or 0.20 (fine), well short of the value 
2NcNf signifying that lattice saturation artifacts have set 
in - indeed for a = 0.186fm trends in all observables look 
smooth out to /ia = 1.0 corresponding to /i = 1.06GeV. 
It is therefore plausible that the observed difference in 
is physical, and due to the differing temperatures of 
the two lattices used. 

Between and jj, d the system resembles the quarky- 
onic matter recently postulated on the basis of large- iVc 
arguments 0] ; namely a state of degenerate matter which 
is also confined, so that excitations above the ground 
state remain color singlet. Because we have used Wil- 
son fermions (with no manifest chiral symmetry) we are 
unfortunately unable at this point to test whether chi- 
ral symmetry is restored, another important aspect of 
the quarkyonic hypothesis; we note however that even 
in a conventional scenario xPT predicts (V'V') oc for 
H > Ho [51] J which is likely to be difficult to distinguish 
from true chiral symmetry restoration in this region, par- 
ticularly if the transition is a crossover. 

The apparent sensitivity of the value oi fiD to small 
changes in temperature is consistent with the very weak 
curvature of the phase boundary between confined and 
deconfined phases postulated in that work, and observed 
in a recent study of QC2D matter using the PNJL 
model [10|. 

An interesting issue is whether QC2D is special in that 
the A^c-quark bound states required by color confinement 
are also favoured by the more general renormalisation 
group argument that 2-body interactions are the only 
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relevant ones close to a Fermi surface [ll|. QC2D is also 
exceptional, of course, because since (qq) is gauge sin- 
glet there is no CSC phase. Whatever the outcome, to 
our mind the study of deconfinement in this hitherto- 
unexplored physical regime promises to be fascinating. 
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